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ABSTRACT 


By examining the relationship of the individual waves in 
an ocean wave record and the power density spectrum (PDS), 
it is demonstrated that the PDS conventionally obtained using 
spectral analysis can be reasonably approximated, except at 
the very low frequency end, by a pseudo-power density spectrum 
constructed from the individual heights and periods contained 
in the record. The very low frequencies seen in the side 
bands of the PDS are examined as products of amplitude modu- 
lation by analogy with simulated wave records and as an 
indication of the presence of wave groups. In addition, 
wave groups are identified and their properties determined 
in 19 selected ocean wave records by integrating the wave 


record using the procedure of Funke and Mansard (1979). 
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iene ROOUCT ION 


The computer era and the Fast Fourier Transform (FFT) 
analysis that is associated with it permits the manipulation 
of wave records in digital form to produce a power density 
spectrum (PDS) with great efficiency. 

The PDS gives a presentation of the distribution of the 
total wave energy contained in successive frequency bands or 
sometimes in period bands. In this way, the analyst and 
interpreter of wave records has a process that describes the 
more or less random individual wave heights and periods 
present in a wave record with a deterministic approach, where 
the peak spectral frequency or period and its relationship to 
the total energy content are the salient features. 

However, the PDS does not fully describe all practical 
characteristics of the wave field. For example, to synthesize 
a realistic wave field in a laboratory channel, the matching 
of amplitude spectra between the prototype and the simulated 
wave field 1s no longer sufficient for the reason that 
different wave fields may have basically the same frequency 
Spectrum (Funke and Mansard, 1979). 

The particular sequence of high and low waves is an 
additional relevant factor. Therefore, attention to grouping 
Of waves and their successive individual properties must be 
included as an additional objective of the analysis of wave 


records. 


eS 





Also, it will be shown that the two evident parameters 
which can be detected by an observer when studying sequences 
of individual waves, their heights and periods (see defini- 
mom tn Section III.A.2), can be further manipulated and 
utilized to generate a so-called pseudo-spectrum. This 
pseudo-spectrum may be considered another descriptor of 
energy in relation to the individual waves that are seen in 
the field. 

The general objectives of this thesis, therefore, are to 
inquire into the Poison wor eorGGaucing a pseudo-frequency 
Spectrum from the height and period measured from the indi- 
vidual waves constituting the wave record, and to analyze 
sequences of heights and periods. 

With regard to the relation between spectrum form and 
mem properties, it is also an objective of the thesis to 
maguire into their detection if possible through considera- 
tion of apparent peaked symmetrical side bands occurring in 
the PDS which may represent the effect of low-frequency 
amplitude modulation. 

Miteelyye te 2S OUL intention in” the course of the study 
presented by Funke and Mansard (1979) to examine wave group 
properties, namely their "degree" of presence in the wave 


Becord, their duration and their variability. 
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It. WAVE DATA 


Zee DATA SOURCE 

The body of this work was accomplished with data from the 
Point Santa Cruz wavericer buoy obtained during the period 
from June 1978 to February 1979. The location, shown in 
Figure 1, is in a water depth of 60 meters. This is deep to 
intermediate water for all ordinary gravity waves (sea and 
swell) as defined by the relative depth. 

The wave buoy was installed as part of the California 
Coastal Engineering Data Network (CEDN) which is sponsored 
jointly by the California Department of Boating and Waterways 
meer tne U.S. Army Corps of Engineers. Additional funding for 
the wave buoy system is supplied by the University of Cali- 
Perma Sea Grant Program with Scripps Institution of Oceano- 
graphy providing overall direction for the actual data 
Seerrection program. 

The data from the buoy received several times a day is 
recorded on a magnetic tape compatible with the IBM 360/67 
computer system. Each record consists of 1024 discrete points 
of instantaneous surface elevation at 1 second interval. 

Thus each wave record has a duration of 17 minutes and 4 
seconds. 

The individual wave records are spectrally analyzed by 


CEDN and the results are published monthly. Tabular data, 


ee, 





an example of which is shown in Figure 2, include: 

Significant wave height 

Total energy 

Percentage of total energy in period bands of 2 seconds. 
Also a visual presentation of the daily energy spectra is 
given so as to provide an indication of the overall daily 
energy pattern for a one month period (Figure 3). This 
presentation enables one to visualize on what days of the 
month the arene content is high, how the energy is parti- 


tioned between the shorter and longer period waves, and how 


broad or narrow is the spectrum. 


B. WAVE RECORD ANALYSIS 

For this study wave records were provided on a magnetic 
tape by the CEDN. Spectral and other analysis of the data 
were performed independently of the CEDN computations. The 
published analyses made by CEDN were used only for purposes 
of comparison with certain of our own analyses. The wave 
records selected for analysis in this thesis were chosen for 
their high total energy and for the dominantly unimodal 


character of their spectra. 


fee LF LLE NUMBER 

For convenience in discussing individual wave records 
and the products of all analyses, file numbers were assigned 
to all records in agreement with those given by CEDN. It 
Should be noted that the figures are arranged in a logical 


order but that they are not necessarily discussed in the 
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following sections in the order in which they appear. Because 
of the large number of printouts contained in the 19 files 
analyzed the decision was made to include only one complete 
file in the thesis for illustration (File Number 228). The 
wave record for this file, which was considered a good repre- 


sentation of a typical wave field, is shown in Figure 4. 


D. CHARACTERISTICS OF WAVE RECORDS ANALYZED 

The significant wave characteristics of the 19 wave 
records analyzed are given in Table I. The wave records are 
classified in the order of the wave type as defined by their 
significant steepness (after Thompson and Reynolds, 1978), 


where steepness is defined as 


27 Hg 
g Tp? 


Sumenalogy with linear theory, H. 1s significant wave height 
and Th is peak spectral period. The defining steepness 


limits adopted were the following: 


Sea 1f/i2 - 1/40 
Young swell 1/40 - 1/100 
Moderate swell Ot — si 250 
Old swell > 7 250 . 
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Iii. INDIVIDUAL WAVE ANALYSIS 


This chapter considers the question of whether the 
Frequency spectrum of a wave record can be approximately 
constructed from the H and T measures of the individual waves 
constituting the record. If this is shown to be feasible, it 
mmeeelt that it will permit a fuller interpretation of the 
frequency spectrum in terms of the waves composing it. 

The approach taken in addressing this question was to 
produce a true spectrum from a given wave record using the 
FFT technique and also a so-called pseudo-spectrum using the 
H and T measures from the same record. 

Bemevarlous reasons, spectra are difficult to compare. 
Accordingly, both spectra were converted into cumulative 
energy distributions and these were then compared. These 
procedures and the principal results obtained are described 


below. 


A. SPECTRAL ANALYSIS 


i True Spectrumt 


1 Although the values are “spectral estimates" of the 
"true spectrum", with uncertainty described by the Chi-squared 
probability density function with eight degrees of freedom, 
the name true spectrum was chosen as contrasted to the "yseudo- 
Spectrum’ derived from the individual H; and Tj. 
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a. Power Density Spectrum (PDS) 

With a basic program that utilizes the Inter- 
national Mathematical and Statistical Library subroutine 
RHARM, the power density was found via the Fast Fourier 
Transform performed on each record of length 1024 points 
with one second interval. 

The data were divided into four "windows" which 
results in spectral estimates distributed according to a 
Chi-squared density with eight degrees of freedom, while 
Maintaining sufficient resolution to detect the low-frequency 
components for periods as long as four minutes. Thus the 
width of each window is 256 seconds, the maximum period that 
can be detected. Consequently, the resolution is 1/256 Hz 
ames 00S9 Hz. 

It is worthwhile to note that for more "windows" 
of shorter duration, say 8 with 128 seconds width, the number 
of frequency components in the PDS decreases. This increases 
the stability of the spectral estimates, since in the fewer 
but broader spectral bands the number of degrees of freedom 
is increased. For comparison, Figures 10 and 11 show PDS 
produced for the same wave record (File Number 228) using 
Four and eight windows, respectively (note that the energy 
density scales differ). 

b. True Energy Histogram 
Also using a class interval of Af = 1/256 Hz a 


histogram of the distribution of energy in percentage of 


Za 





the total was computed (Figure 12) from the PDS (Figure 10). 
This true energy histogram represents the distribution of 
energy centered about the same frequencies that define the 
PDS. In a sense each value represents the area of a "slice" 
resulting from the sum of two pieces with base width Af/2 
extending on each side of the central frequency. 

For comparison, an energy histogram with class 
interval 1/128 Hz produced from the PDS obtained using eight 
Maemeows (Figure 11) is shown in Figure 13. In both figures, 
these histograms are labelled true energy. The values are 
connected by a line as an aid to visual inspection. 

eve cumulative True Energy Distribution 

Barly trials showed that comparison of the true 
energy with our pseudo-energy, not described yet, was better 
done with a cumulative curve, summed from low to high 
frequency. 

The normalized cumulative energy distribution 
computed from the energy histogram in Figure 12 is shown in 
Figure 15 (labelled true spectrum), with ordinates of 
percentage of the total energy. 

A cumulative energy description in energy units 
of length squared was also produced. Figure 14 presents 


this curve also labelled true spectrum. 
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2. Pseudo oreo er ull 


los ie noationeak Ha and Ty 

Using the zero-upcrossing method, the period, 
T;, and the height, Hy, of successive individual waves were 
measured in a wave record. The still water level was 
obtained from the mean of the 1024 instantaneous surface 
elevation digital points (the zero level reference). The 
individual period, Tj, was taken as the time difference 
between two Successive zero-upcross points. The individual 
Memeemenecicit, H;, was taken as the crest to trough distance 
where the crest and trough were represented by the maximum 
positive and minimum negative values, respectively, within 
Beem Pair of Zero-upcross points. 

b. Pseudo-Energy Histogram 

The energy histogram that results from a power 
density spectrum (PDS) is a histogram of energy versus 
frequency, where the energy is interpreted as energy per 
unit area of sea surface per unit weight, with dimensions of 
length squared. This energy is associated with the central 
frequency of the band. 

It is also possible to treat a sequence of waves 
so that individual Hj and Tj are presented in terms of Hs? 
and 1/T;, the same dimensions as above. It is logical that 
the average of the H; in the wave record should represent 
a parameter that is directly related to the energy per unit 


area because those are the waves in that space. 
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Thus if each H,¢ divided by the total number of 
waves, N, in the record lends its contribution to the average 
total energy per unit area, we can use it as an energy para- 
meter for waves within a class interval as well. Then, the 
same class intervals previously used to produce the true 
histogram can produce a pseudo-energy histogram. 

Pisce CimeaPpprGach, the imdividual wave period 
was converted to its reciprocal. This parameter, with the 
dimensions of frequency, was then organized into class groups. 


Hy? 


Then the summation of in each class interval was computed 


and this value was found in percentage of the total sum in 
all classes. With the central frequency taken as reference, 
a plot of the summed values in each class interval yielded a 
histogram. Figures 12 and 13 present histograms, labelled 
pseudo-energy, plotted using class intervals of 1/256 Hz and 
1/128 Hz, respectively. These are superimposed for comparison 
on the true energy histograms derived from the PDS. Both 
sets of curves are produced from the same wave record shown 
in Figure 4 (File Number 228). A comparison of the two 
curves is discussed later. 
c. Cumulative Pseudo-Energy Distribution 

In order to compare more objectively the pseudo- 
energy histogram derived from H; and T; values with the true 
energy histogram derived from the PDS, the histograms were 


transformed into cumulative energy distributions. This was 
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done by cumulating the energy in each class interval (frequency 
band) from low to high frequency in the same manner as with 

the true energy histogram. Figure 15 provides an example of 
the cumulative pseudo-energy distribution in percentage for 
File Number 228, derived from the vseudo-energy histogram 

with a 1/256 Hz class interval (Figure 12). 

With the objective of standardizing the total 
energy in the pseudo-spectrum to the total energy in the true 
Spectrum in energy units of length squared, an attempt was 
made to obtain a conversion parameter. This was accomplished 
by dividing the summation of H; 2 by the number of waves in 
mmemnrecord, N, and by the variance, V, obtained from statis- 


tical analysis, 


ae ae 


Basel 
V 
The interesting result is that a mean value of 6.70 with a 
standard deviation of 0.37 was obtained for the 19 records 
analyzed. These individual values are presented in Table II. 
The reason that the total energy factor appears to be a 
constant and to have the value 6./0 is not clear. 
In view of this apparent constant it was decided 


Sema pply a proportionality factor of 6.7 to convert our 





EO lane wewpal ameter 


individual energy parameter, ra ; 
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Following this, each individual wave was put in ascending 


order by its apparent frequency (the inverse of its period) 


ee 2 
and each individual —_-_— was summed in a cumulative way. 


The results of File Number 228 and of four typical records 
chosen out of 19 records are presented in Figures 14, 32, 42, 
52, and 62 in comparison with the cumulative true energy 
meet htoution. 

eee comparison of Pseudo and True Spectra 

Since the presentation of the cumulative energy 
distribution in energy units depends on the proportionality 
factor chosen, the cumulative curves were compared in 
percentage (see examples in Figures 15, 32A, 42A, 52A and 62A). 

The significant points that result from comparing the 
two cumulative energy curves derived from the 19 wave records 
analyzed are presented below: 

The two curves agree quite closely about the central 
part of the distribution, but differ most toward both ends. 
Thus the slope of the central portion of the pseudo-spectrum 
identifies the true spectrum peak and also characterizes the 
proper spectral width. 

The true spectrum in all of the records analyzed 
identifies a significant amount of energy in frequencies 
lower than the frequency equivalent to the reciprocal of the 
maximum wave period detected by the zero-upcrossing method. 


At this frequency it was found that the differences in 
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percentage of cumulated energy between the true and pseudo- 
spectra among the 19 wave records ranges from 1.75 to 18. 
Thus the lower frequencies of the true energy distribution 
are poorly represented in the pseudo-energy distribution. 

It is also evident that the pseudo-spectrum differs 
substantially from the true spectrum, in the sense that it is 
less smooth. This means that the same quantity of energy is 
described with differently weighted components by the pseudo- 


Seeccrum. 


B. WAVE HEIGHT, PERIOD, AND STEEPNESS 

As already described, the production of the pseudo- 
cumulative spectrum starts with the computation of the 
successive wave heights and periods. It was of some interest 
for the purpose of visually comparing wave heights and periods 
with both the true and the pseudo-spectra produced, and also 
in regard to the examination of wave groups occurring in the 
wave record (described later), to print out in graphical form 
the sequence of Hj and Ty occurring in the record. From our 
19 files we present the sequences of heights and periods for 
Five wave records in Figures 5, 6 (File Number 228), 24, 25 
(File Number 579), 34, 35 (File Number 600), 44, 45 (File 
Number 446), and 54, 55 (File Number 677). 

In addition, the individual wave steepness derived from 


Hj; and Tj was also computed for each wave and the sequence of 
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steepness values was also printed out graphically for each 
wave record (examples are shown in Figures 7, 26, 36, 46, and 
Bo)”. 

The individual wave steepness, defined by analogy with 


the linear wave theory, is given by 


as L 
Ss: = = 


2 
Lj g Ty 





\ 


We are aware that in real water waves the length of an 
individual wave may be different from the length given by 
linear theory. Nevertheless, we reasonably assume that the 
Meeeey ©, Cefined here 1s proportional to the height-to- 
length ratio in real waves. 

Finally, we present in Figures 8, 27, 37, 47, and 57 
graphs of H, versus Tj; these can be interpreted as a 
representation of the individual wave steepness distribution 
in the record. Also plotted on each graph is the significant 


ee SCECepNness, 


27H 
= 2 
g T, 


—" 
(I 


which characterizes the wave record. 
1. Sequences of Individual Height and Period (H,, T,;) 
The following is a summary of the most significant 


points that were found from examination of the sequence 


Seeens Of H; and T;: 
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The sequence of wave heights gives a first indication 
about the presence of wave groups and permits one to count 
the number of waves ina group (see Figures 35, 45, and 55). 
It is interesting to note that the number of successive waves 
that exceeds the significant wave height is normally one to 
two and rarely three to four or more waves. It would be 
interesting to compare these observations with the theory of, 
e.g-., Ewing (1973) on the length of runs of high waves as 
a function of spectral bandwidth, but that is beyond the 
scope of this thesis. 

In the sequence of periods, the spectral peak is 
consistently higher than the mean period, as expected. 
However, in records that can be characterized as "sea", the 
difference is small (Figure 24), and the sequence of periods 
appears to be random; whereas, when looking at records of 
young and moderate swell, such as shown in Figures 44 and 54, 
well-defined wave groups can actually be distinguished by 
their higher values and patterned sequence. 

feeseauence Of Individual Steepness, Sj 


It was found that records classified as "sea" (Table 
I) display more uniform values of steepness (Figure 26), 
whereas steepness values in young and moderate swell have 
higher variability (see Figures 36 and 46; note the different 
vertical scales). Figure 56, as expected with old swell, 


presents a preponderantly constant steepness. All sequence 


Beapns Of S; Characteristically show steepness spikes, 
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3. Individual Height Versus Period 


Each plotted point in these graphs represents the 
steepness S; of an individual wave in the wave record. The 
points are seen, in most of the graphs, to fall within a 
generally well-defined steepness envelope. The envelope is 
associated with large steepness values for sea and young 
swell (Figures 27 and 37) and low values for old swell 
(Figure 57), as expected. In most graphs the width of the 
Gumemope, representing the scatter of the S;, values, is 
relatively narrow (Figure 27) but in several it is broad 
Seagure 66). 

In some graphs the data points show a pattern suggesting 
more than one steepness envelope. In Figure 8, for example, 
the data points associated with the highest and steepest 
waves are contained within the steepness values of 1/40 and 
1/100, whereas a set of lower, less steep waves lies within 
the envelope bounded approximately by the steepness values 
of 1/100 and 1/250. In Figure 57, the higher waves have the 
lower steepness, their envelope being centered about Sj = 1/250, 
whereas, a second set of waves having low heights but high 
steepness appear to occupy an envelope contained between the 
steepness values of about 1/40 to 1/50. 

The significant steepness of the wave record which 
is also shown on each graph is generally representative of 


the waves of greatest height and period in the record, as 
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may be expected. However, it was found that for a few records, 
Significant steepness is not a good representation of the 
waves observed (Figure 66). This situation may be due to the 
presence of a young sea superimposed on an old swell, as 


Figure 63 suggests. 


feet OCORRELATION OF SURFACE ELEVATION 

The application of the autocorrelation function to a wave 
record and the analysis of its periodicity is, as is known, 
another way of finding the main spectral frequency of the 
record (see Figures 9, 28, 38, 48, and 58). However, this 
process can also provide an indication of the presence of 
wave groups. 

As an example, compare the autocorrelation of a sea record 
(Figure 28) with an old swell record (Figure 58). The small 
exponential decay in the swell correlogram is in correspondence 
with the form of the envelope of the wave record (Figure 53), 
and is a distinct "Signature" of the presence of well-defined 
wave-groups, in contrast with the rapid decay observed in the 


sea type wave record. 


D. COMPARISON OF TOTAL ENERGY PARAMETERS 

In the analysis of wave records to obtain the cumulative 
true-energy and cumulative pseudo-energy distributions several 
estimates of the total energy were obtained. Later another 


Function designated SIWEH, but not described yet, was 
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obtained whose mean value is also an indication of the total 
energy. With the objective of comparison, Table III presents 
for the 19 wave records analyzed the total energy computation 
obtained using the following procedures: 

- Statistical variance 

Bpeemcegration of the PDS 


- Mean of the Smooth Instantaneous Wave Energy History or 
SIWEH (Section IV-B) 


pmeumuLative pseudo-energy using a factor of proportionality 
of 6.7 (Section III-A-2(c) ) 


- Total energy presented by CEDN 

From this table it can be said that the first three of 
these methods of computing the variance from a wave record 
present almost equal results, the differences being due to 
the differences in the methods of numerical computation. It 
is of interest that the cumulative pseudo-energy computation 
using a factor of 6.7 gives a value for the total energy with 
mmmaverage error of plus or minus 5.3 percent in relation to 


the variance. 
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IV. WAVE GROUP ANALYSIS 


The question examined here about ‘what process may be 
interpreted as operating in the peak spectral frequency was 
due to the graphs showing the sequence of H; and T; that were 
produced for the identification of individual waves and also 
meemecertain characteristics of the PDS. 

Both sequence graphs show high variability and if the 
sequence of wave height as manifested by wave groups can be 
interpreted as the consequence of amplitude modulation, the 
sequence of periods suggests that some kinds of frequency 
modulation may also be present. We have to say to the reader 
that we did not find conclusive evidence for a simple and 
easily identifiable frequency modulation in our wave records 
examined, but the question is raised because it seems that 
more could be done in describing the wave record through a 
more deterministic approach than has been taken before. Thus 
we are assuming that if the PDS of real wave records show 
peaks on either side of the central peak, this should be due 


to amplitude modulation instead of frequency modulation.@ 


2 The possibility should not be ignored that secondary 
peaks can arise as artifacts of the analysis procedure due 
to "leakage". If so, they appear at intervals related to 
the resolution Af. 
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It was felt that a simple mathematical approach to modu- 
lation phenomena will help in the interpretation of narrow 
Bence rPDs and in probing the nature of wave groups in real 
wave records. Accordingly, this is presented in some detail 


meene next three sections. 


A. MODULATION AND WAVE GROUPS 
1. Amplitude Modulation 
In this type of modulation the amplitude of the 
carrier wave, or carrier, is changed in a proportional way 
to the instantaneous value of the modulating signal. We can 


represent this concept by: 


_—_— = 

- — 
= 
he 


{ ° 
[oo | Plas port wave Carrier 
\ 


Modulating signal 


Wave modulated in amplitude 
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For simplification we assume that the modulating wave 
and the carrier are sine waves. In this way we can produce 


something with this form: 


Ec=E, sin emf t el) 
Sse esa th ai (2) 


mM 


eee sn ent t(3) 





We can manipulate expression [3] in this form: 


Pa) if - + eC 


eS Sin 2 Tey oie 


a) 


(Eo + Ey, Sin 2mf,t) sin 27f.t 


EG Sin 20k Ae apo 27£,t sin 27£,t 


Keeping in mind that 


ses ine CoS eos) Tees Weare, 
2 2 
we get: 
EM 
e = Eo Sin eee + = cos 27 (£. an Ey) 
Em 
- z= cos ate tee) 8 
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Mm 


Meee we define as an index of modulation the value m = a 
c 


wemge ts 


, MEG 
ae eo Slt ene ae ~ —z cos Zanes - fue 


ME ¢, 
- —z~ cos enna - ee 


This is the general expression for a wave train modulated in 
amplitude by a single frequency. 

What can be learned from this is that the spectrum 
of such a modulated wave should present three distinct peaks 


Seeeenie frequencies f,, (f, + f es) lam ty et 


ee Bs 


In Figures 68 (lower graph) and 69, we have a simple 
feeaemetorm and its spectrum, with E. = 100 and f, = 0.1 Hz 
fie= 10 sec) modulated by an E, = 80 with f = 0.016 Hz (T = 
60 sec). The spectrum for this simple case can be interpreted 


ma this way: 


= a 
& } 
= 1 
= ; 
— { 
= ' 
] 
1) 
> ~¥ ‘ 
ee ! ' 
x 
= , 
- i“ ae aac 
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The diagram shows a central peak, 


two side bands equally 


spaced from the central peak, and a low frequency component. 


As far as we can understand, the modulating wave in 


true wave records is not as simple as presented, but can be 


visualized by a complex of frequency components. 


In a visual sense, 


it seems possible that the charac- 


emeatlOn Of the low~frequencies, f,, may be achieved by 


extrapolating to the following schematic concept: 
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In this representation it can be understood that the 
existence of symmetrical side bands can be an indication of 
a well-defined low-frequency modulating signal, and thus the 
presence of symmetrical side bands in the PDS may be expected 
to indicate the occurrence of wave groups. So the wave group 
frequency might be deduced from the frequency separation 
between the side bands and the spectral peak. 

me Erequency or Phase Modulation 

In frequency modulation the instantaneous value of 
the carrier frequency is changed by the modulating signal, 
but the amplitude of the carrier is kept constant. Simple 


Modulation of sinusoidal waves is illustrated as follows: 


Modulator 
Signal 


Carrier 


Wave 
modulated 
In 
frequency 
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Here the instantaneous difference between the carrier 
frequency and the frequency of the modulated signal is 
proportional to the instantaneous amplitude of the modulator 
Signal, but is independent of the frequency of the modulator 
Signal. However, the number of times per second that the 
frequency changes around the central frequency is the frequency 
of the modulator signal. Thus the maximum excursion of the 


Rrequency from the central or peak frequency, Af is: 


es eh 
where: Af = maximum change in frequency 
we > emplitude of the modulating signal 
K = constant of proportionality 


To understand the relationship between side bands and carrier, 


mes consider that the initial carrier can be written as: 


e = E sin @(t) 
where: e = instantaneous amplitude 
B = maximum amplitude 
6(t) = angular value at time (t) 


The instantaneous angular velocity is by definition 


the rate of change of 6(t) with time, 


Becca) es 
OF = oe = 27£5 


Thus we can say that our carrier with constant 
frequency f£, has w, = 2mf, and our O(t) = wot +y, where ¥ 


C 


1s the initial phase angle at t = 0. 
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Here, as the frequency is constant 


a Geiss) _ 
ak dt Vo 


Memem corresponds to the f£, (i.e., w, = 27£,). 


c 


Now, according to the previous definition, a sinu- 
soidal signal modulated in frequency by another sinusoidal 
Signal has an instantaneous value of frequency, expressed by: 


££; = Bs ee eeCoomu tr ., 


where: Af = K Em 


and the instantaneous angular velocity should be (multiplying 


both terms by 27) 


Ws =i ate ne AL COS Wet 


1 c 
Miere: W, = instantaneous angular velocity 
Mewes angular velocity of the carrier 
Wn = angular velocity of the modulating signal 
Af = maximum change in frequency 
. DE eeence) 
Now if i eae 
is 
then 6(t) = wm, (t) dt 
O 
which gives we) a — aye Guat a Sin wt 
Win Tei 


therefore the amplitude of our initial wave when under the 


influence of the modulating signal can be expressed by: 


2 Te 


U) 


e = E sin (wit + 
m 


sin oe 
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Based on this expression, the concept of index of 


modulation in frequency is defined by 





i 7 AF = maximum change in frequency 
a fe frequency of modulating signal 
Secepingeinemind that Af = K E_ , which means that 


Af iS proportional to the amplitude of the modulating signal, 
the index of modulation, assuming a constant Af, is inversely 
MeepertioOnal to the frequency of the modulating signal. 

The side bands originated by frequency modulation, 
fegeeae Of a Simple pair Separated from the £, by the value 
of f, as in amplitude modulation, are separated by Af's that 
are multiples of the modulating signal. The number of extra 


side bands depends on the relation between the frequency of 


the modulating signal and Af. In a sense this depends on 
foe )«=s C'inddex of modulation in frequency). 


It can therefore be shown that the development of the 
wave form with frequency modulated by a sinusoidal wave, as 


expressed, can be presented in the following manner: 
e = aa (me) sin wit 
O 
E pf crea teeta) t sin (uy - w)t) 
ih 
E a oo 20 Sim (ao — 20) t] 
Z 


= 


Chap (me) [sin Cw, ~ nw) t Sin (wo - nw) t] 
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where / (me) are the Bessel functions of order n with 
(n) 


argument Mr. 

The objective for this approach is not only to show 
that a typical frequency modulation has these kinds of 
components from a spectral point of view, but also to take 
Menor the fact that for anm; = 0.6 the lateral bands of 
second and higher order are approximately zero because the 
Bessel function with that argument is almost zero. This means 
that only the fundamental and the first harmonic are important 
mmmeecgard to the spectrum of such a signal, so that modulation 
of amplitude and frequency from the spectral point of view 
look the same. 

3. Modulation in Real Waves 

To exemplify the special equivalence described in the 
previous section and with the objective to better interpret 
real waves, the reader is asked to compare Figures 69 and 73, 


derived from the wave records shown in Figures 68 and 7l, 


respectively (lower graphs). As stated, the PDS in Figure 69 
was produced from amplitude modulation with So WO HZ 
emeaet, = 1/60 Hz having ee 100 and E, = 80, whereas 


the PDS in Figure 73 was produced from frequency modulation 
with Whe 0.6 where the Eas se miZeranic se = M7 GOHZ . 
As can be seen, the two spectra are almost identical. Both 
include side bands centered at the same frequencies about 


the fundamental. 
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iteecanebe learned £Erom this discussion that may be 
applicable to real ocean waves is that side bands should be 
expected in narrow band wave records. The presence of side 
bands is probably due to some kind of amplitude modulation 
rather than frequency modulation because the latter inter- 
pretation implies only a very small change in the wave periods 
(see Figure 72) and a constant amplitude, factors normally not 
expected in ocean wave recoréecs. On the other hand, larger 
variability in the periods can be interpreted as the conse- 
quence of a larger me, and necessarily a larger number of 
components should be present in the spectrum. 

Thus it seems more interesting to confirm whether the 
side bands that can be identified in a wave spectrum may be 
associated with a more or less well-cefined amplitude modu- 


lation which appears in the wave record as wave-groups. 


B. SIWEH AND WAVE GROUPS 

When our research on the side bands of the spectrum had 
reached this point a very interesting report by Funke and 
Mansard (1979) on the subject of synthesizing wave records 
mas brought to our attention. 

For the purpose of isolating wave groups in a wave record, 
these authors note that "a more meaningful description of 
group activity could be obtained by computing the distribution 
of wave energy along the time axis where wave energy would 


be defined as the square of the water surface elevation 
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averaged over a period which is a function of the peak 
frequency." In a sense this operation represents a process 
of smoothing the time series; that is, if we look in the 
frequency domain, what they have done is to perform a filter- 
ing process with a low-pass band filter having a cut-off 
frequency equal to the peak spectral frequency. 

The name that they give to the smoothed wave record, the 
Smooth Instantaneous Wave Energy History (SIWEH), seems to 
us to be well chosen and we will use it here. 

ime oseWEH Analysis 

The Smooth Instantaneous Wave Energy History (SIWEH) 

is produced by operating on the wave record using the Bartlett 
window, which is a triangular weighting function defined in 


mathematical terms by the following: 


SIWEH = E(t) 

L ss 3 
mere: E(t) = ree ae (t + tT) .Q,(t) at 

P ESS 

Pp 
eS - T 
fOr = £ TH p 
ais a -T <T< 

and = Q, (t) A |t|/T,, for -T<t<T,, 


QO everywhere else 
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Where: T duration of the finite wave record 


Th 


ae 
p 


For the beginning and end segments of the record we have: 


peak spectral period. 


I 
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The product of the SIWEH operation on the wave record 
is a time series that reflects the envelope of the record, 
and thereby identifies the larger waves that constitute wave 
groups. Examples of SIWEH are shown in Figures 4, 22, 23, 
33, 43, 53, and 63 (upper graphs). Another comment by these 
authors 1S necessary at this point. In their analysis, the 
authors make the important comment that their work was done 
with “paper waves" which means simulated data. Therefore it 
was considered worthwhile here to analyze our thoughts about 
low frequency amplitude modulation utilizing experiments with 
real wave data. 

Also, some parameters were introduced by Funke and 
Mansard to describe wave groups that had not been evaluated 
using real wave data. We took the opportunity here to produce 
such data for our 19 wave records. These parameters and the 


results obtained are described in the following sections. 
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ee. Eine waverGroup Parameters 


The identification of the low-frequency components of 
our wave records was done by the application of the already 
defined SIWEH. 

This analysis of wave groups is based on the concept 
expressed by Funke and Mansard (1979) that it is important to 
characterize at least four parameters, each derived from the 
SIWEH, in order to fully describe wave group activity. 

The first parameter, the groupiness factor, indicates 
whether there are wave groups present and in what degree. 


This factor is defined by 


GF = 





where: E(t) SIWEH 

In words, the groupiness factor is the standard deviation of 
the SIWEH about its mean and normalized with respect to this 
mean. 

The other three parameters describe the average group 
Guration, the group repetition period, and the variability of 
the group repetition period. 

Ne Gabeulate tne Group duration the autocorrelation 
Function was applied to the SIWEH, and compared with the group 


length measured in the following manner: 


- The mean energy E of SIWEH was computed. 
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- Using E as a threshold the intersections between SIWEH 


and E were determined. 

- The slope of SIWEH at the points of intersection on both 
Sides of each wave group was projected to the zero-energy 
line. 


- These intersections with the zero-energy line define the 
Seeeeent of the group. 


For the group repetition period and its variability 
the SIWEH spectral density function was computed using again 
four “windows", which means a maximum detectable period of 
256 seconds with a resolution of 1/256 Hz and eight degrees 
of freedom. 

The results of these computations for the 19 wave 
records tabulated by wave type are presented in Table IV and 
are discussed in the following sections. 

feeevave Group Properties 

Based on the criterion of significant wave steepness 
the 19 wave records analyzed were divided into groups by 
wave type: sea, young swell, moderate swell, and old swell, 
according to the definitions given earlier. 

For illustrative purposes, one wave record of each 


type was chosen as follows: 


File Number 579 Sea 

File Number 600 Young swell 
File Number 446 Moderate swell 
File Number 677 Old swell 


Other illustrations from File Numbers 228, 178, and 
451 are also presented as examples of particular character- 


istics that will be pointed out later. 
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a. Groupiness Factor 

In Table IV we present the values of the groupi- 
ness factor found. The values vary widely, but show a weak 
relationship to the wave type (see Table I) with sea having 
the highest median value and old swell the lowest value. 

Figure 4 (File Number 228) and Figure 22 (File 
Number 451) represent wave records having high groupiness 
factors. It can be seen that the higher values correspond 
to wave records where a high concentration of energy appears 
after relatively long periods of lower activity close to the 
mean value. It is interesting to note that we did not find 
the high values that are associated with the synthesized waves 
described by Funke and Mansard (1979), for example GF = 1.35. 

b. Duration of Wave Groups 

The duration of wave groups, as expected, acquires 
some meaning only with swell-like records. Here wave groups 
are usually easily identified visually, but examination of 
the SIWEH shows that a typical mean duration of wave groups 
is not easy to define. In sea-like wave records the wave 
groups are not defined and the interpretation of the SIWEH 
has no apparent meaning other than to signify an essentially 
random variation of wave heights. 

The general conclusion from these 19 wave records 
analyzed is that the zero-crossing time lag from the auto- 


correlation graphs (see Figures 16, 30, 40, 50, and 60) does 
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not describe the mean duration of wave groups, but that a 
rough estimative for this parameter can be obtained from the 
eebWinH . 

c. Group Repetition Period and Variability 

The spectrum of SIWEH for swell records gives an 
[iemeation of the periodicity of wave groups (Figures 51 and 
61). The values found and presented in Table IV are condi- 
tioned by our 256 second “window" and the consequent resolu- 
tion of 1/256 Hz. This implies that the peak frequency should 
be a multiple of the fundamental frequency 1/256 Hz such as 
that shown in Figure 20, and this explains the periods 
presented. It may be noted that these periods generally 
coincide with the periods of surf beat commonly seen in 
coastal wave records. 

As discussed above, the presence of wave groups 
may be indicated by the presence of side bands that appear in 
the power density spectrum according to the theory of ampli- 
tude modulation. Thus Figure 59 of File Number 677 shows two 
Side bands present in the PDS that are symmetrically distri- 
buted about the peak spectral period of 15 seconds. This 
low-frequency Signal, also visualized in the envelope of the 
wave record, should have a spectral profile identical to the 
Side bands. The PDS of the SIWEH, Figure 61, confirms this 
characteristic. A similar situation is shown in Figures 10 


emer l/ of File Number 228. 
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We recognize that it is not evident from the PDS 
of the sea-surface elevation which spectral component of the 
Side bands represents the peak frequency such as was produced 
with the simulated waves (Figures 68 and 69). Thus it is 
important in defining the group repetition period to identify 
this low frequency signal, and the application of the SIWEH 
is one process for doing this. 

Thus, as the PDS of the SIWEH of the simulated 
waves (Figure 70) shows the group repetition period of the 
Simulated wave groups, the PDS of the SIWEH of File Number 
677 (Figure 61) presents a group repetition period that agrees 
with the occurrence of wave groups visualized on Figure 53 
(upper graph). 

The variability of the wave group period can be 
seen in the bandwidth of the spectral representation of the 
SIWEH. A broad bandwidth is an indication of large varia- 
bility and a narrow one an indication of small variability. 
In Table IV we categorize the variability only as large or 
small because of the difficulty of establishing some quanti- 
fied parameter. It may be seen that the variability of wave 
group period tends to be large with young swell and small with 
Old swell. The form of the PDS of the SIWEH presented in 
Figures 17, 31, 41, 51, and 61 is instructive when compared 
with the simulated example presented in Figure 70 where 


theoretically the variability 1s Zero. 
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Vv. _CONCLUSZONS 


This thesis presents the results of experiments in the 
analysis of ocean wave records. For this purpose 19 digitized 
records obtained from a wave buoy in 60 meters depth off the 
California coast were analyzed. Two general approaches were 
taken: (1) examination of the interrelationship of the 
properties, H and T, of the individual waves in a wave record, 
and (2) study of the properties of the larger sets of waves 
in a wave record, termed wave groups, and consideration of 
wave groups as a manifestation of amplitude modulation. The 
principal findings are as follows. 

Wave heights and periods, determined using the zero-upcross 
method and treated as a function of H? versus 1/T can be 
Manipulated to produce a cumulative pseudo-energy distribution 
that approximates the cumulative true power density. Thus it 
is possible to reasonably approximate some important proper- 
ties of the PDS by the use of the H and T values obtained 
from a wave record. 

Sequences of H and T for each wave record were produced 
graphically in order to provide different ways of looking at 
the wave records visually. The presence of groups of 
successive large waves is readily seen in the wave height 
sequence plots. Examination of period sequences plots show 


these larger waves to have almost constant and longer periods. 


a 





Plots of H versus T were also produced. Those give an 
indication of how the individual wave steepness is scattered 
about the significant record steepness. 

As a@ manifestation of amplitude modulation, the side 
bands in the PDS of the sea surface elevation give for narrow 
band spectra (swell type waves) an indication of wave group 
Beiivity and periodicity. 

The Smooth Instantaneous Wave Energy History (SIWEH), 
developed by Funke and Mansard (1979) from work with 
synthesized waves, and the power density spectrum of the 
SIWEH were found to provide a good additional means of 
characterizing wave groups in narrow band ocean wave records 
and to give information on the frequency of occurrence, 
Variability, and duration of wave groups. They proved less 
useful for studying wide band spectra (sea waves). 

The groupiness factor derived from the SIWEH was found 
to be a fair descriptor of the degree of development of wave 
Meeuos. tithe application of the autocorrelation function to 
the SIWEH was found to give inconclusive results about the 


duration of wave groups. 
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TABLE II 


TOTAL ENERGY PROPORTIONALITY FACTOR 
BETWEEN PSEUDO AND TRUE SPECTRA 
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FIG 3. EXAMPLE OF CEDN GRAPHICAL SPECTRA PRESENTATION. 
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FIG 20. POWER DENSITY SPECTRUM OF SIWEH IN TABULAR FORMAT (4 windows). FILE NBR 228. 
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FIG 21. POWER DENSITY SPECTRUM OF SURFACE ELEVATION AND SIWER IN TABULAR FORMAT. (8 windows) 
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